N94- 34181 


Gerald L. Morrison, Robert E. DeOtte, Jr., 
Purandar G. Das, and H. Davis Thames 
Texas A&M University 
College Station, Texas 


ECCENTRICITY EFFECTS UPON THE FLOW FIELD INSIDE A WHIRLING ANNULAR SEAL 

rY J 

f- / 

ABSTRACT 

The flow field inside a whirling annular seal operating at a Reynolds number of 24,000 and 
a Taylor number of 6,600 has been measured using a 3-D laser Doppler anemometer system Two 
eccentricity ratios were considered, 0.10 and 0.50. The seal has a diameter of 164 mm, is 37.3 
mm long and has a clearance of 1.27 mm. The rotor was mounted eccentrically on the shaft such 
that the whirl ratio is 1.0 and the rotor follows a circular orbit. The mean axial velocity is not 
uniform around the circumference of the seal; near the inlet a region characterized by high velocity 
is evident on the pressure side of the seal, but the velocity is relatively uniform by the center plane 
of the seal. By the exit, another region of high axial velocity has developed, this time on the 
suction side of the seal. The magnitude and azimuthal distance of the migration increased with 
increasing whirl amplitude (eccentricity). Throughout the seal length the azimuthal mean velocity 
varied inversely with the mean axial velocity. Increasing the whirl amplitude did not increase the 
magnitude of the azimuthal velocity at the seal exit. 

NOMENCLATURE 

_ Tangential mean velocity 
Tune averaged Reynolds stress 
Rotor surface speed 
Distance from seal inlet 
Eccentricity ratio, e/c 
Azimuthal angle 
Fluid absolute viscosity 
Kinematic viscosity 

Fluid density 

Whirl ratio, u p /a> r 
Rotor precession speed 
Rotor speed 

INTRODUCTION 

A seal, by definition, separates regions of high and low pressure fluids. In some cases 
the desire is to prevent leakage, but that is not the general situation. Usually the seals in rotating 
machinery, such as pumps and compressors, are designed to control but not eliminate the leakage 
ow which may cool and enhance operating stability by improving vibration damping. Seal design 
involves balancing considerations of machine efficiency, which is generally higher with lower 
leakage, with operational stability. 

ii / n T y ap P lications - the shaft itself Plays the part of the seal. The seal clearance is 
u ter f ° f 3 millimeter ' Obtaining accurate measurements in these small regions 
is, at best, difficult. As a consequence, the database of velocity and turbulence measurements is 
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very small. Further complicating matters is the turbulent flow found under normal operating 
conditions. The lack of detailed flow field information hinders the development of accurate seal 
performance analyses. 

Research in this area has been mainly devoted to studying leakage and rotordynamic 
coefficients. Though quite important for gaining an insight into the behavior of these seals, they 
alone do not present a complete picture. Velocity measurements have been made in models with 
clearance ratios (seal clearance to seal diameter) much larger than those found in actual 
applications. Turbulence measurements, however, are much scarcer. Most of the data available 
is for seals operating under static conditions and does not provide the knowledge of both the mean 
velocities and the turbulent quantities necessary for numerical modelling of the flow fields. 

Although ideal operating conditions are assumed for seal design and modelling purposes, 
this is often not the case. During machine operation, shafts become misaligned resulting in the seal 
spinning eccentric to the center of rotation thus altering the dynamic forces. The flow in such an 
eccentric seal is different from that in a concentric seal. To provide data applicable to such routine 
situations, the rotors/seals studied were deliberately spun off center. 

One of the important objectives of this study is to help build a database of flow 
measurements for annular seals over a range of operating conditions. Earlier work done at this lab, 
Thames[1992], focussed on flow through an eccentric seal with a 0.5 eccentricity ratio. That study 
represents an extreme case of misalignment. This work compares an annular seal, operating with 
a 0.1 eccentricity ratio, to the 0.5 eccentricity ratio case operating at the same conditions. 

Literature Review 

Nelson and Nguyen[1987] working at Texas A & M University, studied the effects of 
eccentricity on rotordynamic coefficients. They presented a new analysis procedure that solved for 
flow variables in an eccentric annular seal. The governing equations were based on a turbulent 
bulk flow model and Moody’s friction factor equation. Comparisons with experimental data led 
them to conclude that this method was more accurate in predicting the rotordynamic coefficients 
than one which used a finite element analysis procedure. The accuracy of the procedure could be 
improved significantly if a more accurate friction factor formula for annular seals was found. 

Chen and Jackson! 1987] studied the influence of eccentricity and misalignment on the 
performance of high-pressure annular seals. Using existing concentric seal theories they developed 
a method for the calculation of both seal leakage and dynamic coefficients for short seals with large 
eccentricities and/or misalignment of the shaft. The Space Shuttle Main Engine High Pressure Fuel 
Turbopump pump interstage seal was used for calculating the leakage coefficients. From the 
ensuing predictions they concluded that the effect of eccentricity on the leakage rate was higher in 
the case of laminar flow than when the flow was turbulent. They also found that as the 
eccentricity/misalignment ratio (inlet versus outlet) increased the mass flow leakage decreased due 
to an increase in the passage blockage. 

Allaire, et al.[1978], while studying the dynamics of short eccentric plain seals with high 
axial Reynold’s numbers, saw a "Bernoulli" effect in eccentric seals. A bulk flow theory 
developed by Hirsf 1973] was used in the analysis. They found that for a given axial pressure drop 
across the seal, the axial fluid velocity is greatest in the maximum clearance portion of the seal and 
least in the minimum clearance portion of the seal. The pressure distribution around the shaft was 
found to be inversely related to the velocity distribution. This pressure distribution induced a 
hydrodynamic force which tended to center the shaft. At low Reynold’s numbers the pressure 
differences in the seal were small and mainly due to friction effects. At high Reynold s numbers 
the friction effects were small and the pressure was constant around and along the seal. In 
conclusion, they stated that the load capacity of the seal was maximized for a given value of the 
axial Reynold's number, when the Bernoulli effect was the greatest. At very high Reynold’s 
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clearances whereas for the no pre-swirl and backward pre-swirl cases substantial flow development 
occurred for each clearance. Pre-swirl does have a pronounced effect, although it is not 
significantly greater for the statically eccentric than the concentric configuration^ 

Thames [1992] measured the flow field inside a whirling annular seal (e-0.5)usinga 3- 
Laser Doppler Velocimeter with phase averaging. Two axial flow conditions (Re- 12 (XX) and 
R e =24 000) were considered. Velocity and Reynolds stress distributions were presented for both 
flow conditions. The existence of a region of high axial momentum, which at the inlet is on the 
high pressure side of the seal and migrates to the suction side at the exit of the seal, was identified. 


OBTECTIVES 

Most studies of seals have been empirical or numerical. The time averaged Navier-Stokes 
equations were solved by finite difference techniques using a x-e turbulence closure mode . 
Although the numerical solutions predict the flow field, experimental data «s required for 
verification This study was undertaken to measure the flow field through eccentrically whirling 
annular seals. Measurements made were of the instantaneous velocity components in three 
independent directions. From these, the mean flow velocity m the axial radial and tangentia 
directions were calculated along with the turbulent kinetic energy and the phase averaged Reynolds 

stresses. 

FACILITIES AND INSTRUMENTATION 

The experiments were performed in the seal test facility of the Turbomach.nery Laboratoty 
which has been described previously by Morrison, et al.[1991a]. The test section (Figure 1) 
consists of a 50.8 mm diameter overhung shaft which is rotated by a variable speed electric motor. 
The rotor is manufactured of acrylic and is mounted onto the shaft via a bushing. For the present 
study mounting bushings were made so that the inner hole was eccentric with the outer surface 
such that eccentricity ratios of 0.1 and 0.5 were obtained. This mounting results in the lotor 
whirling at the same speed as the shaft, i.e. a whirl ratio of 1.0. The stainless steel stator is 
mounted concentric with the overhung shaft. For this set of experiments, the internal contoured 

plug in the plenum (as used by 
Morrison! 199 1 a]) was not 
installed leaving the face of the 
rotor completely exposed to the 
plenum. 

Water is provided to the 
test section from a 19 m 3 storage 
tank via a centrifugal pump, 
throttling valve, and turbine 
meter. The throttling valve is 
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adjusted to produce the required 
flowrate through the seal. The 
water exiting the test section is 
returned to the storage tank. 

During the tests, the water 
temperature remained between 
ITT and 23'T. 

The velocity field inside 
the annular seal is measured using 

a 3-D Laser Doppler Anemometer 

(LDA) system. This system has been described in detail in Morrison, et al.[1991a] and 


Figure 1. Seal Test Facility. 
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towards the azimuthal location of maximum clearance or accelerate and enter the smaller clearance. 
Since a smaller acceleration is required in the axial direction than in the tangential direction, the 
required pressure gradient is smaller and the fluid accelerates axially. Discussion about the 
migration of the axial momentum to the suction side will be delayed until the azimuthal velocities 

have been presented. 

Velocities 

Figure 3 presents contours of the radial mean velocity. The radial velocities are maximum 
in the inlet plane as the flow adjusts from the large plenum cross sectional area to the substantially 
smaller clearance area. The radial velocities become negligible past Z/L=0.21 as the flow adjusts 
to the small channel clearance. The flow sweeps into the seal from the stator wall resulting in the 
large negative velocities at the inlet plane will reach maximum values close to the region where the 
axial velocities are largest. The whirling of the seal tends to block the flow from entering from 
the plenum region upstream of the rotor face. 

Azimuthal Velocities 

The mean azimuthal velocities are shown in Figure 4. Please note that the rotor surface 
speed is 30.9 m/s. The contour levels were adjusted to show more clearly the fluid velocity. The 
inlet plane contours illustrate that the rotor induces a significant amount of preswirl into the fluid 
by the time it reaches the seal inlet. This is more pronounced for the c=0.5 case. The blockage 
of the incoming fluid due to the large eccentricity results is the higher levels of swirl for this case. 
For the e=0.1 condition, the azimuthal mean velocity has a saddle back type profile, as did the 
axial velocity, which is roughly inversely proportional to the axial mean velocity. That is, when 
the axial velocity is high, the azimuthal velocity is low and when the axial velocity is low the 
azimuthal velocity is high. This inverse proportionality is also seen at the inlet for e= 0.5 with the 
azimuthal velocities being larger on the suction side and the axial velocities larger on the pressure 
side. This may be a simple residence time effect with the slower axial velocity increasing the 
residence time of the fluid which allows the tangential shear produced by the rotor to impart more 
azimuthal velocity to the fluid. For the e=0.5 case, the inverse proportionality of the axial and 
azimuthal velocity components continues to the exit plane where the largest azimuthal velocities 
have migrated to the pressure side and the largest axial velocities are on the suction side. This 
migration of the maximum axial momentum flux from the pressure side to the suction side is 
caused over the first half of the seal (0_< Z/L<_0.5) by the fluid expanding from the pressure side 
with its small clearance to the location of maximum clearance. This is most likely due to an 
azimuthal pressure gradient. However, it is the azimuthal velocity which becomes increasingly 
larger on the pressure side pumping additional fluid from the pressure side to the suction side as 
the flow continues further downstream. The same effect is seen in for the c=0.1 case but to a 
lesser extent since the effects of the whirl are decreased by the decreased eccentricity. The mean 
azimuthal velocity approached 10 m/s (32% of the rotor surface speed) at the exit of the seal with 
very little azimuthal variation for the e=0. 1 case but a much greater azimuthal variation for the 
£=0.5 case. This is larger than the 20% measured by Morrison, et al. [1991a] for a non-whirling 
annular seal However, Morrison, et al. also observed inlet azimuthal velocity values of less than 
10% of the rotor speed compared to the 13% to 30% observed in the present case. The elimination 
of the centerline plug in the plenum resulted in significantly higher preswirl for the present study. 
The feet that the same levels of swirl are present at the seal exit for both the e=0. 1 and 0.5 cases 
show that increased eccentricity (whirl amplitude) does not increase the level of swirl induced by 

the rotor. 
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Axial Velocity Variance 
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CONCLUSIONS 

™ e com P' ex ‘hree dimensionality of the flow through whirling annular seals precludes the 
use of such quasi-steady state models as that of a statically eccentric seal ^1^ 7 
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Figure 4. Azimuthal Mean Velocity Contours 
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